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Abstract 

The expansion of the weak couphng phase of two-dimensional QCD on a sphere is 
constructed. It is demonstrated that the phase transition from the weak to the strong 
couphng phase is induced by instantons. A double scaling limit of the theory at the point 
of the phase transition is constructed, the value of string susceptibility is determined to 
be 'jstr = —1- The possibility of instanton induced large phase transitions in four 
dimensional QCD is explored. 
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1. Introduction 

A longstanding problem has been whether QCD has a representation in terms of string 
theory. Recently, this question has been explored, and a string picture discovered, for QCD 
in two dimensions One uses the fact that the partition function and the Wilson loops 
in two-dimensional QCD can be calculated exactly 0,0. These known answers can then 
be interpreted in terms of sums over maps between two-dimensional manifolds, yielding 
the string interpretation. The hope is that some features of the two-dimensional theory, 
which can be understood in detail due to its exact solvability, have their counterparts in 
the real, four- dimensional world. 

Apart from this, the two-dimensional Yang-Mills theory exhibits other amusing fea- 
tures. Douglas and Kazakov 0] considered QCD2 on a sphere and observed that it has a 
large phase transition as one varies the area of the sphere. Unlike other known large N 
phase transitions []5|, which might be regarded as lattice artifacts, this one occurs in the 
continuum version of the theory. The physical reason for this transition was not completely 
clear. 

The string representation describes, by construction, the strong coupling phase of 
QCD. Remember that in QCD2 observables depend only on the product of the area with 
the coupling so that strong coupling is equivalent to large area. It would be interesting to 
understand what happens in the weak coupling phase. Here one might conjecture that the 
string picture could break down. Indeed there could be a phase transition and the weak 
coupling phase might correspond to a theory of point particles, rather than to a string 
theory. To understand this phase it is necessary to construct the 1/A expansion for weak 
coupling. This we shall do in the next section. We find that, starting from the second term, 
all terms in the 1/A expansion vanish in this phase. There are, however, nonperturbative, 
0{e~^^), corrections to the expansion. These corrections are due to the contribution of 
instantons in two-dimensional QCD. In section 3 we calculate the instanton contribution 
to the free energy. We find that while in the weak coupling phase this contribution is 
exponentially small, it blows up as we approach the phase transition point. The transition 
occurs when the entropy of instantons starts dominating over their Boltzmann weight, 
g-Sinat fpj^g density of instantons then goes from 0, for weak coupling, to cxd, for strong 
couplingS. 

^ Other large A'^ phase transitions, those that occur in matrix models, or in the one-plaquette 
model of QCD can also be attributed to instantons @. 
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Indeed, we can therefore argue that instantons are responsible for confinement in 
two dimensions. At first sight this might appear paradoxical. After all confinement in 
QCD2 is usually regarded as a perturbative effect. In fact in any axial gauge the theory 
is quadratic and the gluon propagator gives rise to a linear potential, which yields an 
area behavior for any non-intersecting Wilson loop by summing Feynman diagrams. Thus 
one would say that the large master field is simply Gaussian fluctuations around the 
trivial, classical vacuum = 0, and certainly does not describe instantons. However we 
note that instanton contributions in QCD2 are not easy to distinguish from perturbative 
Gaussian fluctuations. To discuss instantons at all we must compactify space-time, i.e. 
add the point at infinity. Consider QCD2 on a sphere of area A. Then there are classical 
solutions of the Yang-Mills equations which correspond to magnetic monopoles sitting at 
the center of the sphere, giving rise to constant fields, Ff^^, on the sphere. These have 
quantized charges, F^^, ~ n/gA, are topologically distinct form the perturbative vacuum 
and have action S ~ 'n? /g^A. Therefore in the limit of infinite A these contributions are 
indistinguishable from small perturbations about zero field. To ascertain whether they do 
contribute we should study the theory for finite area, calculate the instanton contribution 
and then take the infinite area limit. As we shall see for small area the instantons do not 
contribute, but for large areas they do, in fact they drive the Douglas-Kazakov transition. 

Analogous large instanton induced phase transitions might occur in QCD in higher 
dimensions. We discuss the possibility of large A^ phase transitions for four-dimensional 
QCD and investigate whether they could invalidate a string representation for short- 
distance physics. We present evidence that such transitions are unlikely. 

Finally, in section 4 we explore the behavior of QCD2 near the phase transition, show 
that there exist a doubling scaling limit near this point, construct this limit explicitly 
and find the corresponding "string equation" (that is, the equation which determines the 
specific heat of the model). 

2. The Weak Coupling Phase 

As we have already mentioned, the partition function of a two-dimensional gauge 
theory can be calculated exactly for any manifold M. of genus G and area A. The result 
can be expressed as a sum over all irreducible representations of the gauge group 0: 




(2.1) 



Z{G, AA, A^) = ^(dimi?) 



R 
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with dimi? and C2{R) being the dimension and the quadratic Casimir of the representation 
R (A is related to ^ by A = g'^N). Since the couphng constant and the area of the manifold 
enter Zji^ only as a single combination XA, we will put A = 1 in what follows. We will 
furthermore restrict ourselves to the case of the gauge group U (N) . 

In this section we will analyze the 1/N expansion of the partition function of two- 
dimensional QCD on the sphere, 



Z{G = 0, A, N) = J](dimi2) V^^2(-^). (2.2) 

R 

If A is kept fixed, the free energy of the theory can be expanded as a power series in 1/N: 

_ „ (2.3) 

g=G 

We will show that for A < Ac = tt'^ 



F{G, A, N) ^ lnZ(G, A,N)=J2 N^'^^ffiA). 

g=G 



(2.4) 



for all g > 2. This means that the expansion in the phase of small areas has only 
two terms. Furthermore, we will evaluate the exponential correction to this trivial 1/N 
expansion. 

The partition function ( |2.2|) can be written explicitly as 

(2.5) 



Zo{A,N)= A2(/i,...,/^)e-^^^('i'-''^). 

Ii>...>In 



We have taken into account that the representations R of U{N) can be labeled by a set 
of integers +00 > /i > Z2 > • • • > ^Af > —00, the dimension of such representation being 
given by the Van der Monde determinant 



dimR^ A(/i,...,/ 



N 



1 

h 

N-l 



I 



1 

In 

N-l 
N 



(2.6) 



i<j 



In terms of {/-,} the Casimir operator C2{h, ■ ■ ■ ,In) has the following form 



C2(/i,...,/Ar) =^n,(n,-2i + iV-M) = ^{N^-l) + J2{l 



(2.7) 



where li = Ui — i + N . Note that (2^) is symmetric with respect to any permutation of 
Ij and vanishes when any two of Ij , Ik are equal. Therefore one can drop the condition 
h > . . . > In, and, introducing 1/A^!, extend the summation to aU unrestricted sets of 
{/i, . . . , In}: 

1 +00 TV / \ 2 



+00 



= ^e-^(^^-i) 5^ A2(xi, . . . , x^.) e-^ 

where we have introduced new summation variables Xi = li — . 

In this form Zq{A,N) can be evaluated with the aid of orthogonal polynomials 0, 
which have been widely applied in the random matrix models. The new effects come from 
the discreteness of Xi, which can assume only integer valuesH. 

We introduce a set of polynomials in x, pn{x\a) = x^ + . . . (where we have explicitly 
indicated the dependence on the parameter a = A/2N) orthogonal with respect to the 
following discrete measure: 

{Pn\Pm) = e~'^''^pn{x\a)pm{x\a) = Sl^hnia). (2.9) 

Adding to any row of the Van der Monde determinant ( |2.6| ) any linear combination of the 
previous rows does not change the determinant. This allows us to represent A(a;) as 

A{xi...xn) = det|x^"^|^^^.^_^ = det\pj-i{xi\a)\^.^^ (2.10) 



Then, we expand the determinants in {^^), and use ( J2.9D to carry out the summation over 
Xj. The result is 



2o(^,iV) = -^e-A("-»nft„(a=A). (2.11) 

j=o 

To calculate hn{a) we shall now derive a set of recursion relations they must satisfy. The 
key idea is that, by construction, the polynomial pn is orthogonal to any polynomial of 
degree less than n. This, in turn, implies the relation 

xpn{x\a) = pn+i{x\a) + Sn{a)pn{x\a) + Rn{a)pn-i{x\a). (2.12) 



If A'^ is even, Xi will be, strictly speaking, half-integers. This does not alter our conclusions. 
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Taking the scalar product of both sides of (|2.12| ) with Pn-i one easily derives Rn{c() = 
hn{(y) /hn-i{ci). The same way, it is not difficult to see that Pn-2 and lower degree polyno- 
mials do not enter the right hand side of (|2.12|) . Also, due to the symmetry of our measure 
with respect to a; — > — 5"^ (a) = 0. 

Therefore, it would suffice to derive the relations between and Rn{a). To 

this effect we differentiate the equation 



(2.13) 



to get 



dhn{(y) 
da 



:^e-"^' 



2pn{x\a) 



dpn{x\a) 
da 



X pMa) 



^e x^p1^{x\a) = -{xpn\xpn) 



[2.1A) 



= - (Pn+l + RnPn-l\Pn+l + RnPn-l) 

= — (hn+l + Rnhn-l) = —hn{Rn+l + R-n) 

We have used the fact that {pn \ ) = because is a polynomial of the order n — 2 
and is orthogonal therefore to Pn- Finally, we obtain our recursion relations 



Rn+i{a) = - (^Rnia) + ^ln/i^(a' 
hn+i{a) = Rn+i{a)hn{a) 

They have to be supplemented by the obvious initial conditions 



(2.15) 



(2.16) 



One can eliminate hn{a) to get a closed form of the relation for Rn{a): 

Rn+i{a) - Rn-iia) = --j- InRnia). 

da 



(2.17) 



When viewed as an evolution equation for functions lni?^(a), a being treated as time, 
( p.lTp is called the Volterra equation . 
Once -Rn(ct) are known. 



hnia) = /io(a) JJi2i(a). 



(2.18) 



1=1 



Now we would like to find Rn{oi = A/2N) for all n. Let us notice that a — in the 
N ^ oo limit. This suggests that the leading term of the expansion can be obtained if 
we approximate the sum in ( p.l6|) by an integral. Remarkably, such approximation gives us 
not merely the first term but the whole expansion, all corrections being exponentially 
small. Indeed, we can elucidate this using the Poisson resummation formula: if 



1 



+00 



Fiy) = — fix)e-^y^dx 

ZTT 



then 



Choosing 



we derive 



+00 +00 



J2 F{n)= Yl /(27rm). (2.19) 



m= — oo 



f{x) = e 2'^^^ 



Let us ignore the exponential correction for a moment, then hQ^\A/2N) = y^2TvN/A, and 
the recursion relations are easy to solve with the result 
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In fact, these correspond to the set of Hermite polynomials Hn{x), orthogonal with 
respect to the continuous measure, as if there were no discrete structure at all: 



2 \/ A7TT1 ' 

e-°^ Hr,{x\a)H^{x\a)dx = 5"^— y . (2.22) 



The discreteness of our true measure ( |2.9[ ) reveals itself in the exponentially small cor- 
rections to ho{a). This suggests that hn{a) may be given by ( p.21|) without any 
corrections. We will see that this is indeed the case provided that the number of the 
polynomial n does not exceed a certain critical value, namely, 

n<ncr = N—. (2.23) 



Technically speaking, the reason is that in the process of iteration of ( |2.15| ) the exponential 
error accumulates and can grow sufficiently large. This leads to a phase transition in the 
structure of the discrete orthogonal polynomials at n = ricr- In fact, when n > ricr, Pn{.x\a) 
differ drastically from Hermite polynomials and ( p.21|) is no longer a good approximation. 
To see how this transition occurs, let us trace the contribution from the first correction to 
hQ^a). If we assume the ansatz 



Nn 



\2NJ a 



(2.24) 



where 



we will have 



Nn 
~A 



A 



+ . . .. 



As a consequence of (|2.17 ) one can derive the recursion relation for c^(A), 

2 



n 



The initial conditions for these coefficients are 

co(A)=0, ci(A) 



(2.25) 



87r2iV2 



To solve for c„(A), it will be convenient to introduce a parameter ^ = '^'^ and to define 
the new unknown functions Gn{£) such that 

2 

Then 



(2.26) 



n L 



(2.27) 



and one can easily see that Gn(^) is a polynomial of the order n — 1: 

^1(0 = 1 , G2{i) = - 1) , G^{i) = 2^2 _ ^ 3 ^ 
We can find Gn{0 ^r an arbitrary n by considering the generating functional 

oo 

F(^,^) = 5^^-G.(0 

n=l 
7 



which, by virtue of (|2.27|) , satisfies the partial differential equation 



OF BF 

-2^ 



dz 

The function F{^, z) is fixed uniquely by the requirement that as z 

F{^,z)=z + Oiz^). 
It is not difficult to demonstrate that it is given by 



(2.28) 



z 

e 1-^ , 



{1-zr 

Then all G'n(^) can be determined in terms of contour integrals 



(2.29) 



dz 1 , 



Q 2ni z 



where the contour Cz encircles the point ^ = in such a way that 2; = 1 is outside of the 
contour. 

Let us cast this integral into a more convenient form. Introducing the new integration 
variable 



t 



l-z 



we can represent it as 



(2.30) 



The contour Ct must encircle t = and pass to the right of t = —1. Now, expanding 
(1 + 1)"^ in powers of 1/t we get an explicit expression 



n-l 



(-2)' 



ni 



k=0 



kl {k + l)]{n - k - ly. 



(2.31) 



In fact, 

where L^^2i(a^) are the so-called generalized Laguerre polynomials, defined as the coeffi- 
cients of the expansion 

J2l^'\x)z- = T^-x^e-^^. (2.32) 



n=0 



(1 - z)P+i 



Combining the pieces we find the answer for Rn'. 



Let us emphasize that the second term in ( |2.33| ) must be much less that the first one, 
otherwise our calculation becomes invalid. Evidently, the first term does dominate for any 



finite fixed n as ^ oo. However, a priori, if ti ~ 0{N), ( 2.33 ) may or may not be 
correct. Indeed, since L^^2j^(x) is a polynomial of order n — 1, it contains terms of order 
[A-n"^ N / A)'^~'^ which at n ~ can win over the small factor e~^'^ ^/"^ thus changing the 
behavior of Rn{A/2N). Still, we need to know Rn{A/2N) for n up to n = — 1 to be 
able to calculate the partition function of the two-dimensional QCD Therefore, we 

need to investigate the asymptotics of (|2.33| ) as n ~ ^ oo in more detail. 

To this effect note that the integral in ( |2.30| ) can be evaluated by the steepest descent 
method. Indeed, if n is of order A^, 

GniO = i ^e-^*W 
Jct 

where is a functional which assumes values of order 1, 

47r2 n , / 1 



It is easy to see that the saddle points of are solutions of the quadratic equation 

47r2 A^ 

Recalling the notation, introduced in ( |2.23| ) we find 



t = t± = -i±i./l--^. (2.34) 

V "-cr 

We see that the position of the saddle points in the complex plane is qualitatively different 
for n < Her and n > Ucr- In the first case it is possible to deform the contour as shown in 
the figurei so that the integral is determined by the contribution of the rightmost saddle 
point, t = 



^ Please note that the contours in Fig.l are directed downwards, because the contour Ct is 
assumed to go counterclockwise around ^ = 0. 
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Fig. 1: Lines of equal level for Re and the integration contour in the two 

cases: (left) n < Ucr and (right) n > ricr- 



T=+ioo 1 

rf^giiV<I."(t+)(zr)2 

27rz 



I g-iV4.(t+) 



Now, after some algebra one finds 

n 2t+ + l 



$"(t 



16n 



iVt2^(t+ + l)^ 



£L /l_ _ < 

V Ucr 



and 



$(t+) 



27r2 27r2 
+ 



1 + Jl 



n In, 
1 In 

Tier 2 TT-cr \ — / ][ — 



(2.35) 



A ~ — 



n 



A A 

The quantity we are interested in is, in fact, the combination entering our formula for Rj 
whereB 



-ll 1 



-1/4 o_2 

^ , ,1- — ) e 



(2.36) 



7(x) = Vl — a; ~ t: 111 



1 + 



2 1- 



4s2 - 1 



s=l 



The (— )" comes from the fact that 1 + -q- < and thus $(t+) has, strictly speaking, an 
imaginary part, equal to mnjN . 
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If X = n/ricr < 1 then always 7(x) > 0. This imphes that 
in ( p.33| ), is indeed exponentiaUy smaU, 0{e~^^^). 
Substituting (|2.36| ) into ( |2.33|) we get 



.361), and hence the correction 



Rr 



Nn 



/ \r7 
+ (-) \/^^cr 

TT^ \ Tier 



1 



-1/4 



-7(^) 



+ .... 



(2.38) 



This asymptotic formula is applicable for < n < rtcr- 

Now it is easy to see how the phase transition in two-dimensional QCD occurs. By 
virtue of ( |2.11| ) the partition function depends only on with < n < A^. So, if 
< ficr, all relevant Rn are given by (|2.38| ) and the system is in its weak coupling phase. 
The condition on the critical point is Ucr = N , which amounts to A^r = tt^. This is 
precisely the value found by Douglas and Kazakov. 

Finally, let us evaluate the free energy of QCD in the weak coupling phase. Using 
(|T1|), (|2T8D and we easily get 

F(A, N) =lnZo(A,Ar) 



^_ A(Ar2_i) + iVln/iof— ) 
24^ ^ "V2iV/ 



N-l 



=1 



i) InRA 

' \2NJ 



(2.39) 



24^ ' 



lnA + 



N 



1=1 



Using (|2.26|) and the contour integral representation ( |2.3C1| ) one can, after some manipula- 
tions with the contour integral, reduce it to 



F{A,N) = -^{N'-1) 



— lnA + 2e-^GNiO- 



Using (|2.36|) we obtain our result 
FiA,N) = 



N 



h InA] H V==- 

2 Vl2 / 24 72^2 



A 



A\-i/4 



-7(4) 



+ . 



(2.40) 



(2.41) 



A brief inspection of this formula shows that the 1/A^ expansion is extremely simple: 
all its terms, except for the first two, vanish!. All non-trivial effects are incorporated 
in the nonperturbative corrections to the free energy, which are invisible within the 1/N 
expansion. Still, it would make sense to have a clear physical picture of these corrections, 
as well as the phase transition. A natural hypothesis is that the 0{e~^^^A^^) term in 
( p.41|) is due to the contribution of instantons to the free energy, while the 1/A^ expansion 
(finite in this case) represents the energy of quantum fluctuations around the vacuum of 
perturbation theory. We shall investigate this possibility below. 



^ The 0{N^) term, A/24, is simply the factor in front of the sum (| 
physical significance. 



It does not have any 
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3. The Contribution of Instantons 



In this section we shall evaluate the contribution of instantons to the partition func- 
tion. We recall Witten's analysis of QCD2 [@], where he generalized the Duistermaat- 
Heckman theorem to the non-Abelian case and showed that the partition function is given 
by a sum over contributions localized at classical solutions of the theory. For finite 
Witten showed that the SU{N) path integral is given by a sum over unstable instantons, 
where each instanton contribution is given by a finite, but non-trivial, perturbative ex- 
pansion. The order of perturbation theory required grows with N, and as A — > 00 one 
gets contributions from all orders, even though the theory is so trivial. We will be able to 
indirectly evaluate and sum this infinite perturbative expansion about the instantons. 

First, let us classify the instantons of the two-dimensional gauge theory on a sphere. 
We are still working with the gauge group U{N). By instantons we mean solutions of the 
classical field equations 



which are not gauge transformations of the trivial solution, A^ 
to a gauge transformation) are simply 



0. All such solutions (up 



A^{x) 



(niA^{x) 

n2Af^{x) 




V 

















(3.1) 



nNA^{x) J 



where A^{x) = v4^(0, (p) is the Dirac monopole potential. 



^e(e. 



0, ^^(e. 



1 — cos G 



G and 6 being the polar (spherical) coordinates on S^: and ni, . . . , nfv are arbitrary inte- 
ge J. The a.«o„ „„ding to .uch io.tanto„ .t.aigMfo.wa.d U, calculat J 



5'inst(^l5 



N 



— / d^xtrF^'^F^^ = Ynl = — Y 
4g2 J pA^ ' A ^ 



(3.2) 



1=1 



^ The condition that rii can assume only integer values follows from the Dirac quantization 
condition for magnetic monopole, which is easy to generalize to our case. 



We keep in mind our convention that A = ^ = 1. 
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These instantons are unstable. For SU{N) there are 2N — 2 negative modes. These 
wiU produce a factor of = (—1)^"^ in the determinant of fluctuations about the 

instanton. Nonetheless the exact partition function of two-dimensional QCD is given by 
the sum over all instanton contributions w(ni, . . . , un)'- 

+00 

J2 w(ni,...,nAr)e-^("i'-''^^) (3.3) 

where w{ni, . . . ^un) are determined a perturbative expansion about the instanton. It 
would be hopeless to evaluate these weights directly, but they can be found from the exact 
answer (|2.8| ), by rewriting the sums over the lengths of the rows of the Young tableaux 
using the Poisson resummation formula (note that the magnetic charges are dual to the 
lengths of the rows). For the gauge group U{N) this resummation has been performed 
explicitly by Minahan and Polychronakos 110] i, who find 



.+00 N N J, , 

«;(ni,...,n^)= / J] J] (yf^ - Aiv'^nl) e'^ l^^-i^^ (3.4) 

i=i i<j=i 

where 

Vij — Vi Vj 5 ^ij — ■ 

Let us calculate the contribution of configuration involving only one single instanton 
of charge n = 1, wi = w{ni=i:n2=...=nN=o), and compare in to the contribution of the 
perturbative vacuum, wq = 'i/j(ni=n2=...=njv=o) in the large limit. From 



/-too 
-00 

/ [[dy^A''{y2,...,yN)e ^^^^""^ [[[4n^ - {y^ - yif] 

J —00 ■ o A o 



i=2 1=2 



We can do the integral over 1/2, • • • , yjv by the saddle point method. In fact as far as these 
variables are concerned we have a ordinary Gaussian matrix model of a — 1 dimensional 
matrix, whose eigenvalues are 2/2, 2/3, • • • , yjv- The saddle point values of these variables are 



A similar resummation formula was used by D'Adda et al [11| to argue, as we shall below 



that the large A'' phase transition is due to topologically non-trivial configurations. 
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of order 1 and have their distribution determined by the 0{N'^) term in the exponential, 
IT 'I2yh ^^^^ density of eigenvalues, p{y) is given by 



p{y) 



11 r 



(3.6) 



7T \ A 4^2 ■ 

The product 

N N 

n [47r^ - {y. - yif] = exp{5^1n[47r2 - {y, - ^0^] } 

i=2 1=2 

would modify the exponential by a 0{N) term, which is negligible compared to 0{N'^) 
term, thus not changing p{y). Therefore, in the large N limit we can just substitute the 
saddle point values of yi, i = 2, . . . , N into ( p.6|) , to obtain 

.+00 N N 



> / dyp{y)ln{47r^-{y-yif: 



(3.7) 



n [47r2 - (y, -2/i)2] > _ =exp 
f4 / Piy) 



wherdi 



/ YldyiA'^{y2,...,yN)e~ 



2A Z^i^2 ''i 



N-2 



(0) 



2AJ 



N-2 

n 



2.n!(-) 



(3.^ 



n=0 n=0 

To integrate in ( |3.5[ ) over yi we notice that the contribution to the integral over yi comes 
from the region yi ~ 1/\/N, due to the presence of e~'^^i/2^_ gg we can replace yi by 
in (13) to get 



wi = Zn-1 



-\-oo 



JV 2 

e 2^2/1 (iy^ 



^JV J dyp{y) In(47r2-y2) 



n=0 



Also, comparing (|3.4|) for ni = . . . = tt-at = to (|3.8|) , we see that 

wo = Zn = YI 

n=C 

Then, the ratio of weights behaves at large as 

«;o (A-l)!VAy 



(3.9) 



^ To evaluate Zn-i we use Hermite orthogonal polynomials ( ^.21| ), ( ^.22| ) for a = N/2A. 
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Evaluating the integral and using the Stirling's formula, 




N(l-\nA) 



we obtain 



(-) 



Wq 



where 7(0;) is precisely the function introduced in (|2.37| ). 

Therefore the contribution of the one instanton configuration, relative to the contri- 
bution of vacuum fiuctuations, is given by 



We see that the perturbative vacuum dominates in the partition function in the weak 
coupling phase. The transition at A = n"^ occurs when the statistical weight of instantons 
grows large enough to make them favorable configurations. If we were to evaluate the 
expectation value of a Wilson loop of area Ai on a sphere of area A we would find two 
different results, depending on whether A > or A < n^. One can say that the large 
master field in these two case is different. The small area phase would give a Wilson 
loop that does not confine as we let A tend to infinity. Thus the perturbative master 
field does not give confinement. The correct, confining, behavior of the Wilson loop on 
the plane is obtained by taking the large area phase and then letting A — 00. We thus 
conclude that confinement is a consequence of the large master field containing instanton 
contributions. 

4. Instanton Induced Phase Transitions. 

We have seen that in evaluating the free energy on a sphere of area A for QCD2 the 
instanton density is zero {i.e. A°°, with A < 1) for A < tt^, and is infinite {i.e. A°°, with 
A > 1) for A > n'^. The jump in the instanton density from zero to infinity induces a 
phase transition. The possibility of such an instanton induced transition is present in any 
large N theory. The instanton density is always suppressed exponentially as — ^ 00 due 

To evaluate the factor in front of the exponential in ( |3.1CI| ) one needs a more subtle technique, 
which would be essentially identical to the one we have developed in the previous section. 




(3.10) 



This result is in agreement^ with (|2.41j ) . 
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to the fact that the instanton action behaves as 5'inst ~ 9^ ~ ciN/\, with A = g^N. 
However, the fluctuations about the instanton can give rise to an entropy factor that can 
overcome the Boltzman weight exp[— S'inst]- The dominant term is the contribution of zero 
modes, which produce a factor of [c2/A]^, where A^^ is the number of zero modes, one 
inverse couphng for each zero mode. These appear from the normahzation of the coUective 
coordinates in the coUective coordinate Jacobian. If ~ '^-c^N, which is the case in QCD 
(both in two and in four dimensions), then the density of instantons is proportional to 
[(c2/A)'^3 exp(— ci/A)]^. This density will blow up, as long as C2C3/C1 > e, for A greater 
than the smallest solution of ci/A = C3 log [02/ A]. 

In our case of QCD2 the number of zero modes of the instanton, i.e. the number of 
generators of U{N) that do not commute with the instanton field, is A^^ — (A^ — 1)^ — 1 = 
2N — 2. Thus C3 = 1. Also ci = 2n'^ and C2 = 47r^e. The above criterion is satisfied, 
namely C2C3/C1 = 2e > e, and therefore a phase transition will take place. Therefore, 
if this was the exact instanton density, the transition would occur at A/rr'^ = 0.746729. 
More generally the instanton density is given by [{c2/Xy^ exp(— ci/A + /(A))]^, where /(A) 
arises from the perturbative expansion about the instanton. In our case we were able to 
evaluate this function exactly, see ( p.41|) , and determine the exact instanton density to be 

00 

piA) = i-r-\ — 1--,) -r e , (4.1) 



nN \ J \ A ^ 

the effect of which is to shift the critical to A = n . 

Can such a phase transition occur in four dimensional QCD? Consider putting the 
theory in a box of volume V. Then for small V we can trust perturbation theory as well as 
the dilute gas treatment of instantons. The finite size of the the box will cutoff the size of 
the instantons, thus they will be governed by weak coupling and have small density. For 
large they will therefore disappear. However, as the size of the box increases the relevant 
coupling will increase and we could induce just such a phase transition. This would be 
very similar in nature to the case of QCD2, which does undergo a phase transition for 
finite volume as the volume is increased. 

Similarly, one might contemplate a phase transition for physical observables. Consider 
the expectation value of a Wilson loop of area A. As A is decreased perturbation theory 
becomes more and more reliable. The contribution of instantons will be cutoff since only 
instantons of the size of the loop and smaller will affect the loop, thus again the instanton 
density will be governed by weak coupling. As the size of the loop is increased the density 
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increases. It could be that for N = oo there is again a phase transition. This would imply, 
among the rest, that the string representation of QCD, which is a representation, 
would break down for certain questions, such as the free energy of QCD in a small box, or 
the behavior of small Wilson loops. The string representation of QCD would then only be 
applicable to the infra-red domain. 

To investigate this possibility let us write the instanton density in the dilute gas 
approximation. This calculation was done a long time ago |T^. The instanton density is 
given as. 



f dp , _ or/0.432133e\ ^-+o(x) 

J LV Hp) J 

where the integral is over the scale size of the instanton and 

1 87r2 11 



N 



(4.2) 



A(p) g^ip)N gm 3 



In(Apvp). (4.3) 



This calculation contains the instanton action, the one-loop determinant and the volume 
of the collective modes. Apy is the renormalization scale in the Pauli-Villars regulation 
scheme. The power of the inverse coupling arises from the fact that the instanton is 
embedded in an SU{2) subgroup of 5'C/(A^) and that there are {N^-l)-{{N-2)'^-l)-l = 
4N — 5 gauge zero modes plus 4 translational zero-modes and one scale zero mode. In the 
ordinary vacuum there is no cutoff on the size of instantons, A(p) grows for large p and this 
approximation breaks down in the infra-red. In the situations we are imagining the finite 
size of the box imposes a cutoff on the size of the relevant instantons, which increases as 
the size of the box, or the Wilson loop is increased. For small boxes or loops the effective 
coupling will never be large, and the instanton density will vanish as ^ oo. But as in 
QCD2 we could contemplate a phase transition as the infrared cutoff is removed and the 
effective coupling grows. 

The question of whether or not instantons survive in the large limit was explored 
by H. Neuberger 00. As he noted, in order to have a large transition we require that 
C2C3/C1 > e. In four dimensional QCD 

ci = 1, C2 = 0.432133e, C3 = 2 ^ C2C3/C1 = 0.864267e < e. 



We thank Herbert Neuberger for bringing this reference, which we had forgotten in an earher 
version of this paper, to our attention. 
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Therefore the above requirement is not satisfied and the maximum instanton density be- 
haves as (0.864267)^ ^ 0, as ^ oo. However it is known that the instanton density 
depends on the method of regularization. Different methods of regularization correspond 
in one loop order to a shift in the renomahzation scale fl^, A ^ A', and thus to a mul- 
tiplicative change in the density of the instantons p(A') = (A/A')^^^^"^p(A). Thus if we use 
MS regularization (Apv/Aj^g) = exp(l/12) [1^] and ci is enhanced by a factor of 1.16507, 
in which case the maximal instanton density behaves as (1.00656)^ ^ cxd as ^ oo. Al- 
ternatively if we use a momentum-space subtraction procedure (in Feynman gauge) then 
(Apv/A]y[OM) = 0.37594, C2 is decreased by a factor of 0.166 and therefore the maximal 
instanton density behaves as (0.14378)"^ — > as — > oo. The regularization that makes 
the instanton density largest is the lattice regularization, for which (Apv/Al) = 42.299. In 
fact Neuberger argued that for this lattice regularization instantons do survive as ^ oo. 
His concern was mostly in the contribution of instantons to the /3-function, and the huge 
jump, using lattice A, suggested a change in physics at the scale where instantons became 
relevant. 

The issue of whether or not a large A^ phase transition occurs for physical observables 
as we vary some physical parameter (the volume of space-time or the size of a Wilson 
loop) is surely independent of the regularization scheme. Therefore it is clear that higher 
order calculations of the instanton density, which are also regularization scheme dependent, 
must be performed to settle the issue. To the extent that we can trust the above, one-loop, 
calculation of the instanton density, using the momentum-space subtraction or Pauli-Villars 
regularization schemes, we would conclude that there is no large N instanton induced phase 
transition in QCD^. 

The above considerations are consistent with the idea that in QCD4 there are no large 
A^ phase transitionJll. However due to the regularization dependence, and the absence 
of higher order calculations, it remains an open issue. There is another reason why one 
might suspect that no large A^ phase transition is likely as one varies the size of a Wilson 
loop. More precisely one could have explored a phase transition in a correlation function 
of local operators, say < 0|^'\E'(a;) 'l'\l/(0)|0 >, as a function of x. For small x instantons 
are suppressed and one could ask whether their contribution to this correlation function 
leads to a phase transition as one increases x. This cannot happen if the large A^ theory is 

There might very well be large A'^ phase transitions in lattice QCD as one approaches the 
continuum limit, but these are lattice artifacts. 
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well-behaved, since it would be inconsistent with Euclidean analyticity of the amplitudes. 
In the case of a Wilson loop the situation is somewhat less clear. We do not know of 
the existence of a spectral representation of the loop which would imply analyticity in the 
radius. However the analogy suggests that a phase transition is unlikely. 

If this is the case a QCD string theory is likely to be applicable both for short and for 
long distances. 



5. The Double Scaling Limit. 

Let us now turn to the problem of the double scaling limit at the point of the phase 
transition. It is a typical situation in matrix models that when there is a phase transition 
all coefficients of the expansion will be singular at A = Ac- If the singularity has a 
special type so that the power of singularity in A — Ac grows linearly with the order of 
expansion, g: 



then double scaling limit [|T6| can be defined by 



N^oo, A^Ac] g-^'^ = N^{A- Acf-'^ fixed. 

Such a limit is known to describe certain versions of two-dimensional quantum gravity, the 
expansion becoming the expansion in powers of the string coupling constant ^fgtr, the 
genus expansion of string theorylil. 

In the case of two-dimensional QCD the triviality of the expansion in the weak 
coupling phase prevents the direct identification of such limit. This might not be the case 
if we approach the point of transition through the strong coupling phase, A > Ac- In 
fact, even in the weak coupling phase there is strong evidence that the double scaling limit 
exists. Indeed, the form of singularity in the leading order in A^, found by Douglas and 
Kazakov [||, 

N^f^{A)strong ' N^fl^{A)^eak ~ N\A - A,)^ (5.1) 

implies that we need to fix 

g-^ = N\.^-A)\ 



Note that the string theory we are contemplating is not the string representation of QCD2, 
for which gsn = t/N, but rather some kind of 2d gravity we obtain at the critical point A = Ac- 
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that is, 7str = — 1. 

We can test this conjecture by looking closely at the free energy ( |2.41D . If A — tt^, 

and 

FiA, N) ~ regular terms , 1 ?r e ^ . 

We see that in the exponential disappears so that it takes the form 

' .iV-V3gi/Vp('--i-). (5.2) 



str 



This is the correct behavior of the nonperturbative terms that one expects. Note the 
characteristic exp[— l/f/str]- In fact, ( |5.2| ) is the small Qstr asymptotics of the complete 
nonperturbative expression that will arise in the double scaling limit. 

Let us now try to find this nonperturbative expression in full. To do this we will 
need to know the behavior of the recursion coefficients i?„ in the immediate vicinity of 
the critical point Ac = tt'^. In the following it will be convenient to introduce the scaling 
variable 

(5.3) 



Keeping Qstr fixed is tantamount to fixing x, and taking ^ oo is the same as considering 
TT-cr oo (since Tier = n'^N/A, A being finite). 

We can extract some information about Rn from our asymptotic formula (|2.38|) : 



TT-cr ^ OO, X fixed. 

We have taken into account that 



n ^ X 



n 2/3 
"'cr Tier 



in the double scaling limits. Since the correction to Rn has an alternating sign, it is 
natural to expect this correction to behave differently for odd and even n. 
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Again, this gives large ricr, and large x asymptotics of Rn, like (|5.2|). 
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To see how behaves at finite x, we consider the foUowing ansatz: 
,2 



7 '"cr 



Tier — ^ OO. 



(5.5) 



The subscripts ± refer to the even and odd n, respectively. 

We wiU now write down a system of ordinary differential equations for the functions 



fi±{^)^ f2±{x), f3±{x). This system follows from the Volterra equation (|2.17| ). Replacing 
a = A/2N = 7r^/2ncr we can rewrite ( p.l7| ) as 

2n2 d 



Rn+l — R 



n— 1 



In Rr 



(5.6) 



Then, considering that 



liH\- 



Tic 



1/3' 



we obtain, say, for odd n. 



Rn+i - Rn-l = nir [ fl+\X 



1/3 
Tier 



1/3 
Tier 



/3 



f2+[x+ — 



+71cr( /3+(x J73 

V ^ ncr 



/3 



(5.7) 



-2nU^f[^{x) - 2ner/^+(x) - 2nli\f^^{x) + ^/{'K^)) + . 



On the other hand, to calculate the right hand side of 
remember that Rn depend on ricr both explicitly and implicitly, through x 



' ^dfcln^-' we must 



dx 



drir 



1 1 X 



n=const nll^ 3 Tic 



The final result of this computation is 

--^lni?^(a) = '^-J-lnRM = ^ + ^nlfU^) 
da Tv^ drier 

+ 2n,,(/^_(x)-7r2/i-(x)/(_(x)) 

+ 2ni (/^_ (x) -\{h- {x) + xf[_ (x)) 

- TT^ (x) {X) + /2- (x)/( _ (X)) + TT Vf- (x) /{ _ (x) 

+ .... 



(5.^ 
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(5.9) 



Equating this to (|5.7|) we get 

f[+{x)+f[_{x)=0, 

/^+(x) + /^_(x)=^(/i_(x)+x/(_(x))+7r2(/i_(x)/^_(x) + /2-(x)/{_(x)) 

Similarly, writing down the Volterra equation for even n, it is easy to see that the 
same set of equations will hold if we replace all fj-{x) by fjj^{x), and vice versa. It follows 



then from the third equation (|5.9| ) that 



1 



f!,+ {x) + (X) =- (X) + Xf[_ (X)) + TT^ (x)/^_ (X) + /2- (x)/(_ (x)) 



1 



6- 



1 



(5.10) 



^^{h+{x)+xf[^{x)) +7r2(/i+(x)/2+(x) +/2+(x)/(+(x)) 
-vrViV(x)/(+(x)-i/r:(x) 
which, together withlli /i_|_(x) = — /i_(x), implies 



(5.11) 



Integrating the second equation in (|5.9|) , we derive 



TT 



f^+ + f,_ = - x+—ft_+C 



(5.12) 



where C is the constant of integration. Substituting this into (|5.11| ) and integrating once 
again, we get an equation for /i(x) = /i_(x): 



f^-4xh-^f! = 3D-37v'Ch. 
The constants C and D can be fixed using the large x asymptotics ( |5.4| ): 

y 



(5.13) 



/i(^) 




X 4e 3^ . 



X — s> +00 



TV- 



(5.14) 



The first equation in ( |5.9| ) can be integrated to get fi+{x) = — /i_(x), since both fi+{x) and 
fi-{x) decay exponentiahy at large x, see ( ^.41) . 
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to be C = D = 0. 

Therefore, obeys the Painleve II equation 



TT 



fi - 4x/i - -ft 



0. 



(5.15) 



Now we are in a position to evaluate the free energy in the double scaling limit. 
However, it is more convenient to consider its second derivative with respect to the area 
(the "specific heat capacity") Fn (A) / dA^ . Since 



Fn{A) = N\iiho(a = + ^(AT-n) lni?^(a 



7V-1 



n=l 



it is easy to see that, as a consequence of ( |2.15| ), ( ^.17| ), 

N-l 



and 



dA 



d^F^iA) 



1 



-Rn + ^ Rn 



n=l 



1 



■Rn{Rn+i + R 



dA^ 4A^2- 

With the aid of our double scaling ansatz (^]^) and using ( |5.12| ) we get 



(5.16) 



dA^ 



nz 



1 + ^(/2-(^) + /2+(x) - ^'/f-(x)) +0( — 

fl{x) + 0(— 



Tier 

1-2^1 



n 

Tier 



9 2/3' 
ZTlcr 



Since at A = Ac one has Ucr = N, we see that the part of the specific heat, singular at the 
transition point, is given by 



d^F^jA) 
dA^ 



4/3 
Tier 



f!{x) + .. 



(5.17) 



sing 



The string equation ( ^.15|) and the shape of the specific heat ( ^.17| ) we have ob- 
tained are identical to the ones encountered in the double scaling limit of the unitary 
one-matrix model. This suggests that the theories defined by the two-dimensional QCD 
at the transition point and by the one-plaquette model 0, |1^ (i.e. one-matrix model) of 
QCD may belong to the same universality class. In fact it is possible to construct a very 
special one-plaquette model, adjusting its potential to make it exactly equivalcTit to the 
two-dimensional QCD on a sphere. We will discuss this point elsewhere. 
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6. Conclusions 

We have investigated the 1/N expansion of the weak coupHng phase of two-dimensional 
QCD on a sphere. We demonstrated that the transition from the weak to the strong 
couphng phase is induced by instantons. We found a double scaling limit of the theory at 
the point of the phase transition, and determined the value of string susceptibility to be 

latr = -1- 

We also explored the possibility of instanton induced large N phase transitions in four 
dimensional QCD and presented a calculation that shows that these are unlikely. There are, 
of course, other possible sources of large N phase transitions, but our result lends support 
to the possibility that a string representation of QCD4 could be valid at all distances. 
This leaves us with the conceptual problem to understand how strings manage to behave 
as particles at short distances? On the other hand it offers us the hope of being able 
to test a string representation of QCD4 by comparing it with the known, perturbatively 
accessible, short distance structure of the theory. 
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